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We study the non-Hermitian Hofstadter dynamics of a quantum particle with biased motion on 
a square lattice in the background of a magnetic field. We show that in quasi-momentum space the 
energy spectrum is an overlap of infinitely many inequivalent fractals. The energy levels in each 
fractal are space-filling curves with Hausdorff dimension 2. The band structure of the spectrum is 
similar to a fractal spider net in contrast to the Hofstadter butterfly for unbiased motion. 


I. INTRODUCTION 

It is well known that the energy spectrum of an electri¬ 
cally charged particle moving in the background of an ex¬ 
ternal magnetic magnetic field on a two-dimensional infi¬ 
nite square lattice has a beautiful fractal structure known 
as the Hofstadter butterfly [T] . The energy bands plotted 
against the magnetic flux reveal a complex pattern which 
resembles a butterfly, hence the name “butterfly”. The 
fractality of the Hofstadter’s butterfly is revealed in the 
fact that small regions of the energy spectrum contain 
a distorted copy of a larger region, thus exhibiting the 
self-similarity property at all scales. This property is a 
characteristic feature of a fractal [2]. 

The Hofstadter spectrum was originally obtained in 
a tight-binding model which describes the motion on a 
square lattice of a (spinless) charged particle hopping 
from one site to one of the nearest sites in the presence 
of an external magnetic field perpendicular to the lattice 
plane. 
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FIG. 1. Biased motion in a magnetic field. 


In this paper we address the problem of a biased quan¬ 
tum motion where the motion is constrained - “biased”- 
by the condition that the hopping of the particle to the 
left is now forbidden, while it can still freely hop to the 
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other nearest sites (right, up and down) as indicated in 
Fig.[T} The particle is electrically charged so that its 
motion is affected by the presence of the magnetic field. 
Its dynamic is described by a truncated non Hermitian 
Hofstadter Hamiltonian. 

This quantum system was introduced in Ref. [3] where 
its relation with biased classical random walks - i.e. ran¬ 
dom walks on a square lattice conditioned to move hor¬ 
izontally only to the right, never to the left [4] - was 
studied. Explicit trace identities relating, on the one 
hand, the generating functions of the algebraic area prob¬ 
ability distribution of biased random walks to, on the 
other hand, the traces of powers of the quantum trun¬ 
cated Hamiltonian, were proposed in the so-called com¬ 
mensurate case where the (normalised) magnetic flux per 
unit cell is a rational number. An exact solution for the 
quantum spectrum was derived. It turns out that the 
energy levels depend on two quantum numbers kx and 
ky. The quasi-momentum ky takes continuous values in 
[—7r,7r] and corresponds to the plane-wave Bloch states 
propagating along the vertical axis while the other quan¬ 
tum number kx takes discrete values in the same interval 
[—TT, tt] and fixes the boundary condition on the horizon¬ 
tal axis. 

Following the approach of Ref. [T] where the fractal 
structure of the Hofstadter band spectrum is revealed 
when it is plotted as a function of the external magnetic 
flux, we are going to show below that the energy bands 
of the biased quantum model also exhibit, non surpris¬ 
ingly, a self-similar structure in the same “energy-flux” 
plane. In addition, we will also argue that energy levels 
for different magnetic fluxes reveal - when they are plot¬ 
ted as a function of ky - unusual fractal patterns (“fractal 
carpets”) in the ’’energy-quasi-momentum” plane. 

The paper is organized as follows: in Sec. |n| we con¬ 
sider the non-Hermitian Hofstadter model and its real¬ 
valued energy spectrum. In Sec. [in| we argue visually 
and numerically that the spectrum possesses a fractal 
structure not only for the energy bands plotted against 
the magnetic flux but also for the energy levels plotted 
against the quasi-momentum ky. An analytical argument 
based on Chebyshev nesting is presented in Sec. |IV[ It is 
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followed by yet another analytical argument for the flat¬ 
tened energy bands in Sec. [V| Complex-valued branches 


of the energy spectrum are briefly discussed in Sec. [Vl 
The last section is devoted to our conclusions 


II. THE MODEL AND ITS SPECTRUM 


with p and q naturally relative prime integers, so that 
their greatest common divisor is equal to unity 

gcd{p, q) = l. (6) 

The spectrum is obtained by solving 

Pp,,{E,ky)=e^'^’^^, (7) 


A. The biased quantum model 

The biased quantum model non-Hermitian Hamilto¬ 
nian is 


Hp = T, + Ty + Tp, ( 1 ) 

where the operators and Ty describe the hopping of 
the particle along (the positive directions of) the axes x 
and respectively. Quantum mechanically, the hopping 
operators act on a state at lattice site (m, n) as [5] 


m,n — •> 

Ty^m,n = 


wher4II 


P = 


$ 


( 2 ) 

( 3 ) 


is the flux ^ piercing an elementary plaquette of the 
square lattice counted in units of the elementary flux 
quantum = hcje. The hopping operators obey the 
commutation relation 


where 


Pp,q{E, ky) — 


E — 2 cos ( ky + 27r-r 


( 8 ) 


is a degree q polynomial of the energy E. Solutions of 
Eqs. 0 and determine the possible energies where 
the continuous variables kx G [—tt, tt] and ky G [—tt, tt] 
play the role of quasi-momenta and the magnetic flux /3 
enters in Eq. (§ via the integers p and q which label 
distinct energy branches. 

The spectrum consists of the q eigenenergies 


Eq r{ky) = 2 COS 


1 

Q 


^arccos 


cos{qky) 


^iqkx 

2 


■ 27rr 


( 9 ) 


where r = 1, 2,..., g' maps different branches of the so¬ 
lutions. It includes, in general, complex valued energies 
which correspond either to formal instabilities or to a 
dissipative motion of the particle. It turns out that the 
spectrum is real for at least certain intervals of the quasi¬ 
momentum ky if and only if kx is such that 

^ (10) 


T^Ty = e-^^^^TyT^ . (4) 

The absence of the T~^ operator in Eq. Q - corre- 
spending to the absence of the hops to the left according 
to Eig.f^- makes the Hamiltonian non-Hermitian be¬ 
cause tJ = for i = x^y and, consequently, ^ 

The normalised magnetic flux [3 enters only via Eqs. 0 
and Thus, the energy spectrum is periodic with re¬ 
spect to the integer shifts, ^ /3 -h 1. By periodicity, we 
take (3 in the interval [0,1]. 


which corresponds to periodic = 1 (antiperiodic 

^iqkx — boundary conditions along the g-site-long 
lattice cell in the x direction. Below, we restrict our¬ 
selves to periodic and antiperiodic boundary conditions 
for which real eigenenergies do exist. 

An interesting feature of the spectrum is that it is in¬ 
dependent of the integer p provided p and q are relative 
prime integers as specified in Eq. 0. 

2. Examples and symmetries of the speetrum 


B. The energy spectrum 

1. General solution 

In Ref. [3] the spectrum of the Hamiltonian 0 was 
derived in the commensurate case where the magnetic 
flux Q is a rational number 


^ Here we use the notations of Ref. In the notations of Ref. [3] 
7 = 277/3 has been used instead. 


In Eig. (top) and (bottom) for g = 7 and 8 respec¬ 
tively, the real energy levels are plotted as a function 
of ky both for periodic/antiperiodic boundary conditions 
(solid/dashed lines). One sees that the q real branches 
labelled by r = 1, 2,... g in Eq. 0 do materialize for 
particular intervals of ky. 

Since the features of Eig. are generic for odd and even 
g, we describe them in details below. 

Eor odd g, the highest/lowest real energy level is a 
continuous curve (solid/dashed) for the whole quasi¬ 
momentum interval ky G [—7r,7r]. The other g — 1 
branches are real only in g distinct ky intervals (here we 
use the fact that the spectrum lies on a circle thus the 
quasi-momenta ky = — tt and tt are identified). These 
g — 1 energy branches are grouped pairwise to form 
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the whole interval ky G [—tt, tt]. The other q — 2 branches 
are real only in q distinct ky intervals. The q — 2 energy 
branches are grouped pairwise to form {q — 2){2 distorted 
ovals in each of these intervals. Thus one has two contin¬ 
uous energy bands and q{q — 2)/2 distorted energy ovals 
(see in Fig. (bottom) the 24 solid ovals for g = 8). 
The antiperiodic case is a bit simpler since the q energy 
branches are grouped pairwise to make g/2 separate en¬ 
ergy ovals in the q bands. Thus one has g^/2 distorted 
energy ovals (see Fig. (bottom) the 32 dashed ovals for 
q = 8). 

The symmetries of the even-g spectrum are simpler: 
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FIG. 2. The energy spectrum Eq for odd g = 7 (top) and 
even g = 8 (bottom) as the functions of the quasi-momentum 
ky for periodic (solid lines) and antiperiodic (dashed lines) 
boundary conditions on the x axis. The g energy branches 
are marked by different colours. 


(g — l)/2 distorted ovals in each of these intervals. Thus, 
both in the periodic and antiperiodic cases, one has one 
continuous energy band and g(g — l)/2 distorted energy 
ovals [see in Fig. (top) the 21 energy dashed/solid ovals 
for g = 7]. 

The odd-g energy spectrum has a number of symme¬ 
tries which are manifest in Fig. i (top): 


(ky) = i-ky)= (ky + 

oGodd \ 


27rn\ 

Q ) 


, (11) 


with n G Z. Indeed, the spectrum is symmetric un¬ 
der both mirroring in the quasi-momentum direction, 
ky fhe discrete shifts ky ^ ky 27rn/g. 

Moreover, for odd g, the periodic and antiperiodic spec¬ 
tra are related to each other by the mirroring trans¬ 
formation in the energy direction E —E and a si¬ 
multaneous shift along the quasi-momentum direction, 
ky ^ ky ^ 7r(2n + l)/g (here we use again that in the 
quasi-momentum space ky and ky + 27r are identified). 

For even g, in the periodic case, both the highest and 
the lowest real energy levels are continuous curves for 


iky) = i-ky) = iky) 


qEeven 




( 12 ) 


with n The spectrum is symmetric with respect to 
mirroring of both energy E —E and quasi-momentum 
ky —ky. The spectrum is also invariant under the dis¬ 
crete shifts of the momenta, ky ^ ky 27rn/g. Contrary 
to g odd, for g even the energy spectra for periodic and 
antiperiodic boundary conditions are not simply related. 

Finally, due to the independence of the spectrum on 
the numerator p of the (normalised) rational magnetic 
fiux 13 in (|^, the spectra in Fig. (top) (bottom) are 
realised, in fact, for six different values of f3 = pjl with 

p = 1,.. .6. 

As clearly seen in Fig. the real spectrum occupies 
the domain [— tt, tt] (g) [—2, 2] with energy bands depend¬ 
ing on the magnetic fiux /3. Following Ref. [3] we plot 
in Fig. I^the energy bands as a function of /3 = 1/g for 
g = 1,2,..., 100. Notice that the band spectrum does 
not carry any information on the actual quasi-momentum 
dependence since it is obtained by a Id projection of the 
spectrum on the vertical (energy) axis. In the next sec¬ 
tion we are going to show that not only the /d-dependence 
of the band spectrum but also the quasi-momentum ky- 
dependence of the energy spectrum do have fractal 
self-similar patterns. 


III. STRUCTURE OF ENERGY LEVELS 
A. Structure of energy bands 

It is well-known that the Hofstadter band spectrun0 
E G [—4,4] reveals a fascinating fractal structure - the 
Hofstadter butterfiy - when it is plotted on the horizontal 
axis against the normalised magnetic fiux /3 G [0,1] on the 
vertical axis. Likewise, in the biased quantum model, a 
fractal structure should also materialize when the energy 
bands are plotted against 


^ Notice that in the biased quantum model the real-valued band 
spectrum lies in the narrower interval E G [—3, 3]. 
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FIG. 3. The band energy spectrum (the horizontal axis) with 
respect to the minimal value of the normalised magnetic flux 
/3 = 1/g G [0,1] with q = 1 , 2 ,..., 100 (the vertical axis) for 
periodic boundary conditions. 



In Fig. 1^ the band spectrum is plotted against P with 
all values p/q E [0,1] (and not only the minimal values 
1/q as in Fig. for g = 1,..., 70, both for periodic and 
antiperiodic boundary conditions. Here q plays a role of 
a “resolution” parameter which determines the thinness 
of the chosen grid of possible P values. 



FIG. 4. The band energy spectrum (the horizontal axis) 
with respect to the value of the normalised magnetic flux 
P — p/q ^ [0,1] with g = 1, 2,..., 70 and p = 1,..., g (the ver¬ 
tical axis) for the periodic (the red lines) and antiperiodic (the 
blue lines) boundary conditions. The darkness/lightness of 
the lines visualise high/low values of the density of states ( |13| ). 


of states at a given energy 


1 dky 

27r dE 


(13) 


such that the higher/lower density of states, the 
darker/lighter the band. As an illustration in Fig Jm (top) 
and (bottom), we display the density of states ( p^ and 
the corresponding energy levels for g = 2,4,8, respec¬ 
tively. 





FIG. 5. (top): the density of states ( p^ as a function of the 
energy E for g = 2 (the blue solid line), g = 4 (the orange 
long-dashed line) and g = 8 (the violet short-dashed line) 
levels; (bottom): the energy levels as a function of the quasi¬ 
momentum ky. 

At a fixed resolution g the (normalised) magnetic flux 
runs over the values 


q q q 


(14) 


where P = 0 corresponds to p = g. The energy levels are 
defined by Eqs. 0 and ^ that is to say 


'[[[E-2 cos {ky + 27r/3r)] = . 


(15) 


r=l 


In Fig.j^the variations of lightness in each energy band 
correspond to the variations of the (normalised) density 


Evidently, some of the instances of p and g in Eq. (14) 
correspond to mutually reducible integers so that their 
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greatest common divisor gcd(p, g) 7 ^ 1. In this case one 
can use the property of the energy polynomial ^ which 
appears in the LHS of Eq. ([T^ 


[E — 2 cos {ky + 271Pr)] 


r=l 

where 




/5=f 


Qq{E, ky) = 2 [Tg {E/2) - cos{qky)] , (17) 

and Tq{x) is a Chebyshev polynomial of the first kincQ 
Substituting Eq. (16) into Eq. (15), we get for the en¬ 
ergy levels 


which can be simplified as 

Qg^{E,ky)=e^i^’^% 
where the integer qp is 

Q 


(18) 


(19) 


( 20 ) 


m 


gcd(p, q) ■ 

We come to some trivial conclusions for the spectrum 
the case of the rational (normalised) magnetic flux ^ 

1 .) if p and q are mutual prime integers as in (|^, then 
the energy spectrum is determined by 


Qg{E,ky)=e^‘>'^^ 


( 21 ) 


because qi = q. 


2 .) if p and q are not mutual primes [gcd(p, g) ^ 1 ], 
then the energy spectrum is determined by 


Q<7p ^v) ~~ ^ 


\qkx 


( 22 ) 


with qp given in Eq. (20). 


3.) in the special case p = q^ where the magnetic field 
is absent, /3 = 0 , we get the simple g-independent 
solution 




o'^kx 


cos{ky) + — 


(23) 


which is Eq. for g = 1 and thus r = 1 . 


Eor p = 1,2,...g—1 the integer qp runs over all possible 
divisors of the natural number g. Therefore the nontrivial 
energy levels at resolution g are determined by Eq. (22) 
where the g^’s are taken from the set of all non equal 
divisors of g. 


^ We use the standard notation Tn for the Chebyshev polynomial 
of the first kind determined by the equation Tn(cos cp) = cos(n(p). 


Eor example, the scan of non zero possible normalised 
magnetic fluxes with resolution g = 6 includes the se¬ 
ries P = 1/6,..., 5/6. According to our discussion, the 
mutually prime numbers in (3 = 1/6, 5/6 share the com¬ 
mon spectrum corresponding to /3 = 1/6, while /3 = 2/6 
and /3 = 4/6 correspond to the spectrum of ^ = 1/3 
and P = 3/6 corresponds to the spectrum of P = 1/2. 
Therefore, the magnetic fluxes at resolution g = 6 should 
include the solutions at [3 = 1/2,1/3,1/6 [plus, of course, 
the trivial solution ([^]. In other words, the solutions 
are given by Eq. with g = 1, 2,3, 6. 

As yet another example, consider the magnetic fluxes 
with resolution g = 16, given hy ^ = 1/16,..., 15/16, 
which can be divided into four different groups. The first 
group with the mutually prime numbers, P = (2n+l)/16, 
n = 0,1,..., 7, has a spectrum corresponding to the 
minimal representative p = 1/16. The second group, 
given hy p = 2{2n + 1)/16 with n = 0,... 3, corre¬ 
sponds to the mutually prime numbers P = (2n + l)/8 
with the minimal representative P = 1/8. The third 
group is P = 4/16,12/16 which has the same spec¬ 
trum as P = 1/4. Einally, the last group includes only 
P = 8/16 = 1/2. Thus, the magnetic fluxes at resolu¬ 
tion g = 16 should include the solutions corresponding 
to P = 1/2,1/4,1/8,1/16 and (23). In other words, the 
solutions are given by Eq. in which g = 1,2,4, 8,16 
runs over the set of all divisors of 16. 

Both examples illustrate that the energy levels with 
resolution g are given by the the energy solutions (§ 
in which g^ runs over all possible divisors of the original 
resolution factor g [including, of course, the trivial energy 
level Eq. (23)]. 

The energy bands are plotted in Eig. they do have 
a characteristic pattern which no longer resembles a but¬ 
terfly but rather a spider net. 


B. Fractal energy levels in quasi-momentum space 

1. Self-similarity 

To address the fractal structure of the eigenenergies of 
the non Hermitian Hofstadter Hamiltonian 0 plotted as 
a function of the quasi-momentum ky one has to consider 
curves made of all the disjoints energy ovals lying in the 
planar domain [— 7 r, 7 r] (g) [—3,3]. 

Let us focus on the series g = where the base m in¬ 
cludes any integer 2,3,5,6,7,10,... which is not a mem¬ 
ber of another series with a smaller base (i.e. m ^ 
for any integers mi > 1 and ni > 1). Eor example, Eig.l^ 
(top) for g = 7 is the starting plot of the series q = 7^ 
with base m = 7 whereas Eig. (bottom) for g = 8 be¬ 
longs to the m = 2 series g = 2^4^8^16^32^ 

.... We are going to argue that for a given base m in the 
limit n ^ oc the curve made of the reunion of all the 
disjoint ovals in each of the g = members of the se¬ 
ries constitute a fractal area-filling curve with Hausdorff 
dimension 2 , as it is the case for well-known area-filling 
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curves such as, for example, the Peano, Wunderlich or 
Hilbert curves [2]. 

The curves for the series q = m." with m = 2,3,5 are 
shown in Fig. Note that they are not per se strictly 



2 

1 
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FIG. 6. The energy spectra § are shown as functions of the 
quasi-momentum ky at (normalised) magnetic fields — Ijq 
for q = mT' series with (a) m = 2, n = 1,..., 6; (b) m = 3, 
n = 1,..., 4; (c) m = 5, n = 1,..., 3. 

self similar because of the distortion of the ovals near 
the upper and lower energies E = ±2. Leaving aside 
these spurious boundary effects, the overall pattern is 


those of a fractal which remains unchanged by a more 
an more accurate zooming. In the sequel we concentrate 
on the central part of the energy plot, say E G [—1,1], 
where boundary effects become negligible and the self¬ 
similar structure of the oval patterns is manifest. As 
an example, the self-similarity of the energy levels in 
the q = 2^ series is explicitly illustrated in Fig. It is 



FIG. 7. Zooming into the fractal structure of the energy lev¬ 
els of the non-Hermitian Hofstadter model for (normalised) 
magnetic helds ^ —\/q with q — 2^ and increasing n. 

clear that in the central region (far from the distorted 
boundaries) the q = 2^^+^ energy levels are self-similar to 
a zoomed version of the q = 2^ levels, this for any power 
n. This self-similarity depends on the base number m of 
the q = mE series. Thus, in the limit n ^ oo, we have 
not one but rather an infinite number of fractals labeled 
by the base m that one considers. 


2. Energy levels as spaee-filling eurves 

There are several ways to determine the fractal struc¬ 
ture of a 2d curve. One way is to consider a curve which, 
by successive iterations, would have a diverging length 
yet the area in which it is contained remains finite (typ¬ 
ical examples of such curves are the Koch and Peano 
curves EI). The scaling of the length of the curve with 
the length’s resolution should give us the Hausdorff - 
fractal - dimension of the curve. Let us first address 
this question numerically and then, in the next section 
Section [TVl consider some more precise analytical argu¬ 
ments. 

As we have already seen, at resolution q the number 
of energy ovals is q{q — 2)/2 for an even q and q{q — 
l)/2 for an odd q. As a simple example, let us again 
consider the q = 7 and q = 8 curves of Figs. Since 
the energy pattern repeats itself horizontally q times, it 
is sufficient to consider the part of the curve made of the 
ovals contained in a given ky interval (a given column). 
Numerical estimates ioi q = 2'^^q = 8'^^q = 7E and large 
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n (typically, the maximal power is n = 24 for base m = 2) 
give for the perimeter length and the area enclosed 
^coi respectively 


^coi ^ 7^34771 _ 


^col ^ 


7.74880 


(24) 


As for the total perimeter length and the total area 
^tot entire curve made of all the disjoint ovals, the 

scaling^ are respectively 

£tot _ 7,34771 X q , - 7.74880 . (25) 


The scaling (25) - a diverging length for a curve which 


nevertheless encloses a finite area - is typical of a fractal 
area-filling curve. The Hausdorff dimension D of the 
curve and the resolution of the curve’s details 1/q - so 
that the bigger the resolution g, the smaller the scale at 
which one looks at the curve (as already mentioned, 1/q 
corresponds to the minimal possible nonzero value of the 
normalised magnetic flux /3) - should be related by 


= 7( - 


1-D 




(26) 


A comparison of Eq. (25) with Eq. (26) gives a Hausdorff 


dimension D = 2 which does correspond to an area-filling 
curve. 

Another way to show that the Hausdorff dimension 
is 2 consists in rephrasing the scaling argument above 
in terms of iterative patterns which are clearly visible 
in Eig. for base m = 2: each time the resolution is 
doubled, i.e. 1/g is divided by 2, the number of ovals is 
multiplied by 4. This property is again characteristic of 
a fractal with Hausdorff dimension D = In4/In 2 = 2. 
In general for base m, iterating means dividing 1/g by 
m with a number of ovals multiplied by m^, so again an 
Hausdorff dimension D = Inm^/lnm = 2. 


IV. CHEBYSHEV NESTING AND FRACTALS 


The energy spectrum of the non-Hermitian Hamilto¬ 
nian ^ is determined by Eqs. ( pT| ), ( [^ or (22) which 


involve Chebyshev polynomials of the first kind. The 
fractal structure of the q = series can surprisingly be 
understood in terms of the so-called “Chebyshev nesting” 
which is a composition identity for Chebyshev polynomi¬ 
als 


Tm{Tn{x))=Tmn{x). (27) 


Note that the area enclosed by the ovals is contained in a rect¬ 
angle of area 27r x 4 ~ 25.1327 whereas the ’’available” surface 
inside the rectangle, i.e. after removing the gaps between the 
bands, is [3] in the large q limit 27r x 8/3 = 16.7552. 


We can use three different approaches to discuss the 


relation of Chebyshev nesting (27) to the fractal self¬ 
similarity of the energy curves. 

Firstly, let us consider the horizontal scaling of the en¬ 
ergy solutions of Eq. (21). As we have already noticed in 
Fig. 13 and Fig. the undistorted self-similar pattern of 
the energy levels is manifest in the vicinity of the E = D 
horizontal axis. The same self-similarity should in partic¬ 
ular be seen in the scaling behaviour of the zero-energy 
solutions which occur for a set of specific quasi-momenta 
ky. In other words, if, for the series q = mE ^ the energy 
curves exhibit a self-similarity as the power n increases, 
then one might expect that the discrete quasi-momenta 
ky for which the energy lines cross the horizontal axis will 


also exhibit a self-similar scaling. According to Eqs. ( |11 
and (21) for periodic boundary conditions - where in¬ 


deed energy lines do cross the horizontal axis for both 
even and odd q - these zero-energy quasi-momenta are 
given by 



(0) = COS 

{rti^ky) + - . 

(28) 

Using the Chebyshev nesting ([2^ as well as 


T,(0)=cos^ = | +1 

for odd g; 

for even g, even g/2; 
for even g, odd g/2; 

(29) 

and 




Tq{El) — +1 , 

T,(-l) = (-l)^ 

(30) 

one arrives at 




T„.(0) 

= T„(T„(...T„(0)))=/+, 

^ > 

(31) 


n 

^ 2 times 



where 


= 

J m 


1 ± (—1)"^ — / 1 even/odd m, 
0 otherwise. 


(32) 


The important consequence of (27) is that, according 


to Eq. (31), T^n(O) does not depend on n for n ^ 2. 
Then, in the periodic case at hand, the zero-energy quasi¬ 
momenta solutions of Eq. (28) are 


^ 0) = ^ T [3 _ (_!)« + 121] , (33) 

where I G Z labels the possible independent solutions 
contained in the interval [—7r,7r]. 

According to Eq. (33), for any fixed two powers n and 


n p with p G Z, the zero-energy quasi-momenta are 
related by 

= 0) = = 0), (34) 


This implies that the zero energy quasi-momenta for q = 
can be obtained by a rescaling of the q = mE quasi¬ 
momenta by a factor 1/m. In other words, the width of 
the ovals diminishes by 1/m with each step n ^ n + 1. 


















Secondly, let us consider the vertical scaling of the en¬ 
ergy ovals for periodic boundary conditions. This scaling 
can be determined in a similar way using a simple argu¬ 
ment based on the Taylor expansion of Chebyshev poly¬ 
nomials. Indeed, let us determine from Eqs. 0 and 
(21) the spectrum in the region close to the E = 0 axis 
where the energy spectrum is undistorted. For a given 
q = mE ^ the Chebyshev polynomial in Eq. 0 can be 
expanded in power series as 

{x) = gm,n{m^x) + 0(a;^) , 


(35) 


where 


gm,n{y) = fm+fm f t-i + (- 1 )"/^ 


y ■ 


(3g) 


de- 


The first two terms in the small-x expansion 
pend solely on mEx. The next-order term in Eq. (|35) 


does not share the same scaling, but in the small ener¬ 
gies region (or, equivalently, small x), it can be neglected 
as well as the higher-order terms. These facts play an 
important role for the fractal properties of the energy 
spectrum: indeed keeping only these lowest-order terms, 
Qm^n is invariant with respect to the shift n ^ n + 2 


9m,n-\-2{y^ — 9m,n{y) 


(37) 


One can then distinguish two sectors according to the 
parity of n 


gZ^^{y) = gm,n{y) 


9 m {y) = 9mAy) 


nGeven 


nGodd 


= A + fmy-^-fy\ ( 38 ) 

= /++2/sin^-^y2^(39) 


Thirdly, there is a bit more elegant and general way to 
show the fractal self-similarity of the spectrum both for 
periodic and antiperiodic boundary conditions by using 
yet another property of the Chebyshev polynomials 


cosx, even m, 

lim [(-l)™"T„.(m-"a;)] = <^ (43) 

smx, odd m. 


For practical purposes it is convenient to use an ap¬ 
proximate formula based on (EF 


T^n (m '^x) -A (m '^x) = cos x 


(44) 


+ (-!)-"/-sinx. 


where is given in Eq. (32). The approximation made 


in (44) is valid in the interval x G [—7r,7r] with an error 
10 “^or m = 2, n = 4. The accuracy of this approxima¬ 
tion is increasing rapidly as either m or n or both become 
larger (for example for m = n = 4 the error is less than 
10-^). 

For large n one can combine, both for periodic and an¬ 
tiperiodic boundary conditions, Eqs. ( p!7[ ), ( [21] ) and (44) 
into 


even m: 


cos 


mX E 
2 

m^E 


= COS {mEky) ± \ , 


odd m: (—sin = cos (m^/c^) ± |, 


(45) 


which are again valid in the energy domain mE\E\ tt. 
One can check that the solutions of Eqs. (45) reproduce 
the lowest energy levels. 


The symmetries of Eqs. (45) clearly imply the fractal 
self-similarity of the energy levels 


where = even (odd) for even (odd) n, respectively. 
For a given 9m independent on n. 

Clearly, for q = mE with periodic boundary conditions, 
the energy close to the horizontal axis is determined by 

5m = cos{m^ky) + i . (40) 


= ^_l)mp^-p^im,n) {jn^ky) , (46) 


both for even and odd m and p G N. These equations 
are consistent with the scalings (34) and (41). 

Equation (46) implies that the energy spectrum for 




q = m 


can be obtained by 


Due to the independence on n of the left hand side of 
Eq. (40) once the parity In is fixed, the energy solutions 
for q = and q = mE are related by 


1.) rescaling (i.e., squeezing) the q = mE solution by 
the factor 1/m both in the energy E and quasi¬ 
momentum ky coordinates, and 


^(m,n+ 2 p) ^ ^ (4^^ 


where p G N is an arbitrary natural number. The self¬ 
similar scaling property (41) is consistent with the zero 
energy quasi-momenta scaling (|M|). 


The explicit solutions of Eqs .^38 j), (|39|), ( |4Q| ), 

—2 cos{m'^ky) 


£j(rn,n) _ ^ 


± 


1+2 cosijnX'ky) 


(- 1 ) 




even m , 


odd m , 


(42) 


are low-energy solutions valid in the region mE\E\ tt. 
They highlight the self-similarity properties ( |M| ) and ( jdT] ) 
of the energy spectrum as n increases. 


2.) periodically copying (i.e., extending) the squeezed 
solution m times along the energy E axis and m 
times along the quasi-momentum ky axis. 

Note that these properties cannot be readily seen from 
the original Eq. They highlight three important 

properties of the q = mE series which were already nu¬ 
merically illustrated above. As n ^ n + 1: 

A) the total area An enclosed by the energy ovals in 
a fixed area of the {ky^E) plane remains constant 
since the ovals become m-times smaller in each of 
the two {ky and E) directions while their number 
increases by a factor so that An+i = An] 
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B) the total lengths and of the energy oval 
perimeters in a fixed column and a fixed row respec¬ 
tively remain constant since they become m-times 
smaller while their numbers in each column and in 
each row increase by a factor m so that 

and £--1 = 

C) the total length of the energy oval perimeters 

in a fixed area in the {ky^E) plane increases by a 
factor m so that due to the second 

property as well as because the number of columns 
and ovals increases by the same factor m. 


When n ^ oc these three properties imply that the 
total length of the energy curves (i.e., of the oval perime¬ 
ters) in a fixed area of the {ky^E) plane diverges as 
(X m'^ while the area A enclosed by these ovals re¬ 
mains constant. Thus, a union of the energy ovals in¬ 
deed constitute an area-filling curve with Hausdorff di¬ 
mension D = 2. In agreement with Eq. (26), the size 
On oc llq = 1 /mE of any individual oval along any direc¬ 
tion decreases with an increasing resolution, g ^ oc, so 


that one gets x a~^ = a^E^ with the usual D = 2 
Hausdorff dimension for space-filling curves. 

Notice that the band structure of the energy levels has 
been obtained by projecting the /c^-dependent energy 
levels on the vertical energy axis and plotting it against 
the the magnetic flux [3. Since the quasi-momentum de¬ 
pendent energy levels are self-similar, we conclude that 
for a fixed base m the energy band spectrum has also to 
exhibit a fractal self-similarity for the different powers n 
of the normalised magnetic flux P = 1 /mE , as illustrated 
in Fig. 


V. FLATTENED ENERGY SPECTRA 


where 


f{x) 


[ 0 , 27 r] 


TT — |x — 7r| = arccos(cosx). (49) 


The flattened spectra (47) are shown in Fig. [^a), (b) 
and (c), respectively for the q = mE series with m = 2, 
3 and 5 - they correspond to the unflattened spectra of 
Fig.[^ Clearly, the arccos-operator in Eq. ( [47| ) “unfolds” 
the spectrum by flattening the energy levels close to E = 
±2, so that no distortion is present anymore. It follows 
that the fractal self-similarity of the spectrum is manifest 
in the whole energy domain. 

For example in Fig. |^a), each of the q = 2^ energy 
ovals intersects with two smaller q = 2^+^ energy ovals 
for n = 1, ... 5. Analogously in Fig. [^b), each of the 
q = 3^ energy ovals intersects with four and includes one 
smaller q = 3^^+^ energy ovals for n = 1,... 3. And in 
Fig. I^c) a similar self-similarity pattern can easily be 
observed for the g = 5^^ series. 


B. Self-similarity for the flattened bands 

Besides these visual observations, one can make a more 
precise statement on the self-similarity. In Fig. we plot 
the flattened band spectra for g = mE with m = 2, 3, 5. 
These plots share certain features with Fig. [^but with 
the difference that in Fig. [^the vertical axis shows the 
power n for a given base m instead of 1/g for a set of g’s. 

Consider first in Fig.j^a) the q = 2'^ series. The spec¬ 
trum at the lowest level n = 0 has - in terms of the 
flattened energy variable @ - a single band of width 
6820 = 27r/3 and a gap of width ir/S. In order to get the 
spectrum at the next level n = 1, it is sufficient to 


A. Flattening 


(i) rescale the original n = 0 spectrum, shrinking it by 
a factor 1/2 so that f G [0,7r] ^ f G [0,7r/2] 


The fractal self-similarity of the real-valued energy 
spectrum can be addressed in yet another way. As 
stressed above, the energy ovals are distorted near both 
edges of the spectrum E ^ ±2 (as illustrated in Fig. [^. 
On the other hand in the center F ~ 0, the distortion 
is practically absent and the self-similarity is quite pro¬ 
nounced. It turns out that in order to remove the bound¬ 
ary distortion effects from the entire energy domain, it 
is sufficient to consider, instead of the energy itself, the 
energy-related quantity 


Sq^ri^y) — arccos 


1 

2 




(47) 


in the interval £ G [0,7r]. 

Using ^ for the g energy branches r = 1,..., g of the 


original spectrum, one can rewrite Eq. (47) 


^q,ri^y)~f 


1 f 

- arccos 

Q V 


cos{qky 


piqkx 


■ 2iTr 


(48) 


(ii) duplicate the resulting spectrum to cover again the 
whole energy range, £ G [0,7r/2] U [7r/2,7r] ^ f G 
[0,7r] 

(hi) and finally invert the duplicate around its center 
£ = 37r/4. 

Obviously, the total width of the two resulting bands 
at level n = 1 remains equal to the width of the single 
original n = 0 band, = ^^ 2 ° • 

Recursively, the band spectrum at level 
n + 1 is obtained by applying the same rescal- 
ing/duplicating/inverting procedure on the spectrum at 
level n: 

(i) rescale the spectrum by a (l/2)-shrinking 

(ii) duplicate the resulting spectrum 

(iii) invert the duplicate around its center £ = 37r/4. 
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FIG. 8. The flattened energy spectra ( [47| vs. the quasi- 
momentum ky for q = series with (a) m = 2, n = 1,..., 6, 
(b) m = 3, n = 1,..., 4, and (c) m = 5, n = 1,..., 3. These 
plots correspond the original energy spectra shown in Fig. 


In Fig. I^a) this procedure is explicitly illustrated up to 
level n = 5. Note that some pairs of adjacent smaller 
bands merge in single bands with a width twice larger. 
Note also that the inversion is not trivial only for the first 
iteration n = 0 ^ n = 1 since from n = 1 onward the 
spectrum becomes symetric around its center. 

As a result, for the q = 2^ series, at level n > 1 there 
are {2^ + 2)/2 bands with individual widths necessarily 


summing to 

< 5 ^ 2 *°* = Y- ( 50 ) 

By construction, this band spectrum is evidently self¬ 
similar, and it becomes a fractal in the limit n ^ oo. 

Likewise, consider in Fig. ib) the flattened energy 
bands for the series g = 3^^ up to level n = 4. Start¬ 
ing from level n the iterative procedure consists in 

(i) rescaling the spectrum by a (l/3)-shrinking; 

(ii) triplicating the resulting spectrum; 


(iii) and finally inverting the second copy around its 
center 5 = 7r/2 

with pairs of adjacent bands again merging into single 
bands of width twice larger. As a result, for the q = 3^ 
series, at level n > 1 there are (3^^ + l)/2 bands which 
again have a total width = 27r/3, as in Eq. (50). 

Similarly to the q = 2'^ series, in the limit n ^ oo the 
q = 3^ band spectrum becomes a fractal. 

A similar iterative procedure applies to the band spec¬ 
trum for the g = 5^^ series as illustrated in Fig. IHc) up 
to level n = 3. The bands at level n + 1 are constructed 
from the bands at level n by a 1/5-shrinking, 5-plicating 
and inversion of every even copy around their centers 
£ = 37r/10 and £ = Ttt/IO respectively. Again adja¬ 
cent bands merge into single bands twice wider so that 
at level n > 1 there are (5^^ + l)/2 bands with a total 
width S£ln^ = 27r/3. Clearly, in the limit n ^ oo the 
g = 5^ spectrum becomes a fractal. 

The universality of the iterative self-similar procedure 
is now evident. For the g = series the band spectrum 
at level n +1 can be constructed from the band spectrum 
at level n by 

(i) (1/m)-shrinking. 


(ii) m-plicating. 


(iii) inverting every even copy around their centers £ = 
(2k — l)7r/(2m) with /c = 2,4,_ 


As a result at level n > 1 there are (m^ + l)/2 (for odd 
m) or (m^ + 2)/2 (for even m) bands with a total width 
d£^n = 27r/3. This implies that the bands always occupy 
2/3 of the available flattened energy range [0,7r]. The 
universality of the procedure is highlighted by the fact 
that all the band spectra are constructed from the same 
initial n = 0 level band spectrum. 

For every given base m each successive iteration n 
n + 1 has the bands shrinking to smaller bands. Thus, in 
the limit n ^ oc, one obtains a point-like spectrum con¬ 
sisting of infinitely many bands of infinitely small width. 
The structure of the band spectrum is similar to a Can¬ 
tor set with a difference, which is, however, essential: 
contrary to the Cantor set, the total band width is the 
nonzero base-independent quantity 


lim 

n^oc 


27T 


( 51 ) 



















11 



(a) 



(b) 



(c) 


FIG. 9. The blue horizontal lines: The arccosine (47) of the 
energy levels § for the series q = mP' at hxed power n and 
for (a) m = 2, (b) m = 3 and (c) m = 5. The shadowed 
regions and the arrows show the self-similarity features. 


Thus, similarly to the original unflattened band spec¬ 
trum, the flattened band spectrum has a fractal self¬ 
similar structure when plotted against the normalised 
magnetic flux P = p/q with p = 1,..., with a reso¬ 
lution q = mP for base m determined by the increas¬ 
ing level number n. The self-similarity pattern depends 
on m but the total band width ( [bl] ) is m-independent. 
The flattened band spectrum has a Hausdorff dimension 
D = 1 which follows from the 1/m-shrinking m-plicating 
iterative procedure. 


Finally, we stress that, as far as the unflattened spec¬ 
trum E(ky) in eq. ^ is concerned, distortion effects tend 
to disappear in the n ^ oc limit: as n increases the en¬ 
ergy “ovals” distortion effects migrate from the center of 
the spectrum E 0 towards its edges, E ±2. In other 
words, in the limit of an infinitely large q = m'^ resolu¬ 
tion, regions with distorted ovals eventually disappear at 
the edges of the energy domain. This feature is clearly 
visible in the plots of Fig. which are worth comparing 
with their corresponding flattened versions in Fig. In 
Ref. [3], the unflattened q ^ oo spectrum was shown to 
have a total band width saturating to 8/3. Compared 
to the total unflattened energy range, E G [—2,2], the 
bands therefore occupy in this limit 2/3 of the whole en¬ 
ergy interval. Equation (51), not surprisingly, leads to 
the same conclusion for the flattened band spectrum. 


VI. COMPLEX ENERGY BRANCHES 


So far we have discussed the real energy branches of the 
non-Hermitian Hamiltonian 0 - However, the full spec¬ 
trum contains energy branches with a non zero imaginary 
part. Since these energies correspond either to some for¬ 
mal instabilities or to a dissipative motion of the parti¬ 
cles, they are not particularly relevant for a large system 
at thermodynamical equilibrium. Still, for the sake of 
completeness, we find appropriate to display the entire 
spectrum with its complex energy modes. 

In Fig. f^the g = 8 spectrum is displayed in a 3d (/c^, 
ReE, ImEy^plot. The real branches shown in Fig. are 
recognisable in the horizontal {ky, ReE, ImE = 0) plane. 
They build ordered chains of distorted ovals. Similarly, 
the complex branches build ordered chains of distorted 
ovals in the {ky^ ReE = 0, ImE) plane. These complex 
(“vertical” in the 3d plot) ovals connect the real (“hori¬ 
zontal”) ovals in a 3d chain-like manner. 

Energy spectra for generic q values, both for peri¬ 
odic and antiperiodic boundary conditions, are similar to 
Eig. with horizontal real ovals always connected by 
vertical complex ovals along the ky coordinate. This ob¬ 
servation gives strong evidence to the fact that the entire 
energy spectrum - with the imaginary branches included 
- is also a fractal, similarly to the real energy spectrum. 
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FIG. 10. The real and imaginary parts of the energy spec¬ 
trum for q = 8 plotted against the quasi-momentum ky. The 
horizontal blue lines represent the real spectrum (ImE = 0) 
while the vertical red-orange lines correspond to the complex 
spectrum with a nonzero imaginary part (ImE 7 ^ 0). A pro¬ 
jection of the energy spectrum on the bottom of the plot is 
also presented for convenience (grey color). 


VII. CONCLUSIONS 

We have studied the fractal properties of the energy 
spectrum of a charged particle with a biased motion on 
a two-dimensional square lattice in the background of a 
magnetic field. This biased quantum model is such that 
the particle is not allowed to hop in one of the directions 
(say, to the left) while it can freely move in all other 
directions. Its dynamics is described by a non-Hermitian 
Hamiltonian of a Hofstadter type. 

We have shown that the energy spectrum, which de¬ 
pends on the normalised magnetic flux - the ratio of the 


magnetic flux threading an elementary plaquette to the 
magnetic flux quantum - possesses a nested multi-fractal 
structure. 

The energy bands, plotted against the rational mag¬ 
netic flux, exhibit a fractal pattern shown in Fig.[^ Con¬ 
trary to the Hofstadter’s butterfly, the band spectrum is 
rather similar to a spider net. 

The energy levels, plotted against the quasi¬ 
momentum, exhibit even more curious fractal patterns 
as illustrated in Fig. We have shown, both numeri¬ 
cally and analytically, that the real energy spectrum is 
an overlap of infinitely many inequivalent fractals which 
we call “fractal energy carpets”. The energy levels in 
each fractal are space-filling curves with Hausdorff di¬ 
mensions 2. 

In a more rigorous approach, we have shown that for 
a given base m and resolution q = with n G N, the 
energy curves for the normalised magnetic flux (3 = 1/ mE 
(m = 2, 3, 5, 6, 7,10,... is required to not be equal to a 
natural power of a natural number) have a self-similar 
geometric structure. This self-similarity is observed up 
to certain finite scales, both in quasi-momentum and in 
energy, constrained by the finiteness of n. In the limit 
of an infinitely fine resolution, n ^ oc, the structure of 
the energy level becomes a fractal with various fractal 
carpets depending on the base m considered. 

We have also shown that the real branches of the en¬ 
ergy spectrum are connected by complex branches form¬ 
ing chain-like structures in the 3d space of the energies 
plotted against the quasi-momentum. 

The fractal properties of the energy levels are visu¬ 
alised in the supplementary video material [6] . 
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